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Abstract 

We present an application that automatically writes the Helas (HELicity Ampli- 
tude Subroutines) library corresponding to the Feynman rules of any in quantum 
field theory Lagrangian. The code is written in Python and takes the Universal 
FeynRules Output (Ufo) as an input. From this input it produces the complete 
set of routines, wave-functions and amplitudes, that are needed for the computa- 
tion of Feynman diagrams at leading as well as at higher orders. The represen- 
tation is language independent and currently it can output routines in Fortran, 
C++, and Python. A few sample applications implemented in the MadGraph 5 
framework are presented. 
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Operating system: Linux/Mac/Windows 
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Classification: 

4.4 Feynman diagrams, 

11.6 Phenomenological and Empirical Models and Theories 
Nature of problem: 

An efficient numerical evaluation of a squared matrix element can be done with the 
help of the helicity routines implemented in the HEL AS library [1] . This static library 
contains a limited number of helicity functions and is therefore not always able to pro- 
vide the needed routine in presence of an arbitrary interaction. This program provides 
a way to automatically create the corresponding routines for any given model. 

Solution method: 

Aloha takes the Feynman rules associated to the vertex obtained from the model 
information (in the Ufo format [2]), and multiply it by the different wavefunctions 
or propagators. As a result the analytical expression of the helicity routines is ob- 
tained. Subsequently, this expression is automatically written in the requested lan- 
guage (Python, Fortran or C++) 

Restrictions: 

The allowed fields are currently spin 0, 1/2, 1 and 2, and the propagators of these 
particles are canonical. 

Unusual features: 
None 

Running time: 

A few seconds for the SM and the MSSM, and up to few minutes for model with spin 
2 particles. 
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1. Introduction and motivation 

We have now entered the era of the Large Hadron Collider (LHC) and mea- 
surements based on the Standard Model (SM) predictions have already begun. 
As data continues to be analysed, deviations from the SM are likely to occur and 
the challenge will be understanding which Beyond the Standard Model theories 
(BSM) are responsible for such deviations. Extracting information from the LHC 
data requires the ability to accurately simulate events based on the theoretical 
models. While it is possible for researchers to manually create event simulators 
for the SM and a few popular BSM theories, such an approach would severely 
restrict the range of models that could be tested at the LHC. 

To alleviate the potential bottleneck in simulating BSM theories, extensive 
efforts are being made to automate the process of going from the Lagrangian 
of a new theory, to simulating events that can be then compared with data 
[HEl El II]- In this paper we present one important piece of that process, i.e., the 
ability to automatically create helicity amplitude routines necessary for calculat- 
ing the matrix element associated with any new model. Such an ability will be 
directly useful to many groups writing custom generators, however, at present the 
greatest impact is on the very recently released MadGraph 5 suite [5]. Prior 
to this work, MadGraph was restricted to generating amplitudes for models 
that had the same Lorentz structures as the Standard Model. Adding new par- 
ticles and new interactions, from the MSSM [6] to higher spin particles such as 
gravitons [7] or gravitinos [8], required identifying, writing and debugging new 
helicity amplitudes associated with the model. This severely restricted the range 
of models that could be studied, as the number of physicists able to implement 
new models is very limited. For example, technicolor, chiral perturbation theory, 
and chromomagnetic moments all require extensions to the HELAS library. The 
ability to dynamically generate the required helicity routines eliminates the bot- 
tleneck, giving all physicists the ability to easily simulate any new model within 
the MadGraph 5 framework. For instance, all models presently available in the 
FeynRules framework [TO] can be used via the UFO format [LI] in Mad- 
Graph 5. 

This paper is organized as follows. In Sec. [2] we give a simple example with the 
aim of explaining the idea of calculating Feynman diagrams routines for vertices 
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and propagators as implemented in Helas. This also allows us to establish a 
dictionary of terms that will be used in the following sections. In Sec. [3] we 
describe in more detail how the ALOHA program works. In Sec. [4] we describe 
the validation procedure and collect our conclusions in Sec. [5j 

2. Helicity amplitudes 

Helicity amplitudes methods [121 US EEU EES! EE CEO EEE] are a convenient and 
effective way to evaluate the squared matrix element of any process. As the name 
implies, helicity amplitude methods work at the amplitude level, in contrast to 
the trace techniques based on completeness relations which work on squared 
amplitudes. This has two important advantages. First, the complexity in the 
calculation grows linearly with the number of diagrams instead of quadratically. 
Second, it is possible to "factorize" diagrams, such that if a particular substruc- 
ture shows up in several diagrams, it only needs to be calculated once, significantly 
increasing the speed of the calculation (see for example [5]). 

The Helas [T9J [20] library has been a particularly successful implementa- 
tion of an helicity amplitude method for tree-level processes. Amplitudes are 
generated by initializing a set of external wavefunctions using their helicity and 
momenta. These wavefunctions are combined based on the particle interactions 
in the Lagrangian to determine the wavefunctions of the internal lines (propa- 
gators). Once all of the wavefunctions are determined, they are combined to 
calculate the complex number corresponding to the amplitude for the diagram. 
These amplitudes can then be added and squared to give the required result. The 
set of routines needed to calculate amplitudes for the Standard Model at tree- 
level was released in the original Helas package [19] and successfully employed 
for many phenomenological studies thereafter. 

In addition to having a complete and highly optimized set of routines, the 
structure of the Helas calls makes it particularly easy to write (and read) the 
code for any Feynman diagram: one line of code for each line of the associated 
Feynman diagram, and a final call that returns the amplitude. For example, the 
diagram for W + W~ — > ti via a Z boson, shown in Fig. [TJ is written in six lines. 
First, there are calls to functions associated with the external particle^] 

call VXXXXX(P_WM, W.MASS, WM.HEL, -1, WM) 
call VXXXXX(P_WP, W_MASS, WP.HEL, -1, WP) 
call 0XXXXX(P_T, T.MASS, T.HEL, +1, F0) 
call IXXXXX(P_TBn T_MASS, TB.HEL, -1, FI) . 

From the four resulting wavefunctions (WM, WP, F0, FI) one can combine either 
the two fermions or the two bosons in order to build the propagator for the Z 
boson. Choosing to connect the fermions, the code reads: 



1 The notation follows the Helas convention |19j . 
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W+ i 




Figure 1: W + W — > tt via Z exchange 

call J3XXXX(FI, FO, GAU, GZU, Z.MASS, Z.WIDTH, J3) . 

Finally, the Z boson can be combined with the external bosons to evaluate the 
amplitude: 

call VVVXXX( WP, WM, J3, GW, AMP). 

The example above demonstrates that Helas consists of three categories of 
routines: 

External particle routines: Each external leg of a Feynman diagram corresponds 
to a wavefunction computed for a given momentum and a given helicity. 
For example, the first four routines used in the W + W~ — > tt calculation 
are of this class. The only input expected for this category of routines is 
the type of the particle (i.e. vector, scalar, tensor or incoming/outgoing 
fermion), is independent of the interactions in the model, and depends only 
on the Lorentz representation of the particles present in the theory. 

Wavefunction/ Off-shell routines: When one leg of a vertex has no associated 
wavefunction, a routine is used to generate the internal particle (i.e., prop- 
agator) wavefunction. This is exemplified by the Z boson call above. These 
routines return both the internal wavefunction and the momentum of the 
propagator. 

Amplitude/ On- shell routines: When the wavefunctions for all the legs at a vertex 
are known, they can be used to compute the matrix-element of the diagram 
by the associated amplitude routines. 

To simulate new models, relatively few external particles are required, as 
they depend only on the Lorentz representation of the particle. However, the 
wavefunction and amplitude routines depend explicitly on the interactions of the 
theory, making it impractical to anticipate and create a comprehensive library of 
routines for all interesting BSM theories. 
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3. Automation of the method 



Originally, the Helas package contained external particle routines only for 
classes of particles present in the Standard Model (particles with spin 0, 1 and 
1/2). Recently the library has been expanded to include spin-2 particles [2"Uj I2"T]. 
spin-3/2 particles [HI E2], as well as the Higgs effective theory [23]. This set of 
external particle routines is sufficient for calculating most renormalizable BSM 
theories of interest. 

The challenge that remains is writing the required wavefunction and ampli- 
tude routines required for any BSM theory, renormalizable or not. Although the 
rules for writing such a routine are known, it can be very time consuming to ver- 
ify that all of the sign conventions and normalizations have been properly imple- 
mented. This provides a severe limitation to HELAS-based event generators. To 
fill this need, we develop Aloha: the Automatic Language-independent Output of 
Helicity Amplitudes, which is capable of automatically generating Helas routines 
in Fortran, C++ and Python. 

3.1. Computation method 

Aloha is written in Python, the same language used by the Ufo code [H] 
to ensure maximum flexibility and compatibility. Moreover, the input model file 
used by Aloha is the one provided by the Ufo. 

In order to generate the necessary Helas routines, Aloha needs to perform 
specific analytic computations. A package to deal with these mandatory compu- 



tations has been developed and its details can be found in |Appendix B 

We begin by describing the relevant information existing in the Ufo code, 
which is the starting point of Aloha. All information concerning the model, 
i. e. particles, interactions and parameters, is obtained from the Ufo. Aloha con- 
nects particles, creating Lorentz structures to output the wavefunctions and am- 
plitudes using the syntax defined by the original Helas package. 

To illustrate how Aloha works, consider the e + e _ 7 vertex produced by Ufo: 

FFV = Lorentz (name = 'FFV' , 

spins = [2, 2, 3 ] , 
structure = J Gamma(3,2, 1) ' ) 

The spins attribute contains the list of spins of each particle (present in the 
vertex) in the 2*S+1 scheme^] The structure attribute collects the analytical 
expression for the Lorentz structure in the vertex. The arguments of the Lorentz 
object refer to the particle associated with each index of the object. The list of 



2 In the case of a fcrmion, the first particle is considered as an incoming fermion and the 
second as an outgoing fermion. For interactions containing more than two fermions, the con- 
vention is that the vertex should be defined for a given fermion flow. The outgoing fermion 
follows directly its associated incoming fcrmion. 
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Table jlj Elementary 


Lorentz structures 
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Right chiral projector: (— 6 


ProjP(a,b) 




Left chiral projector i}-^) b 


ProjM(a,b) 




Sigma matrices: (a' J " u )ab 


Sigma(/i,z/,a 


,b) 


Scalar wavefunction: 0at 


Scalar(N) 




Spinor wavefunction: ijj^ 


Spinor(a,N) 




Vector wavefunction: e„ 


Vector(//,N) 




Spin2 wavefunction: T^' 17 


Spin2(/i,z/,N) 





Table 1: The greek indices stands for lorentz indices while the lower case latin indices stands 
for spin indices. The indices N stands for the label of the particle. 

available objects and the Python convention for different indices are explained in 
Table [TJ The UFO convention requires: 

• positive indices to be linked to the particle numbers; 

• all Lorentz indices to be listed before spin indices; 

• all repeated indices to be summed with the appropriate matrix; 

• all momenta to be incoming. Note this is different for the Helas/ Aloha con- 
vention. 

From a given Ufo model structure, Aloha creates the analytical expression 
linked to the associated Helas routine by contracting the expression with a set 
of wavefunctions. To illustrate how it works we use the structure shown in the 
FFV python object introduced above: Gamma(3,2, 1). 

The associated amplitude routine, which was previously named IOV by the 
HELAS convention, is now: 

Spinor(2) * Gamma (3, 2, 1) * Spinor(l) * Vector (3) 
which correspond to the analytical expression: 

FFV_0(V>i, ifa, e 3fl ,g)=g* u 2 ^ Vl t^. (1) 
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The analytical module of Aloha expands Eq. ([I]) component by component. The 
result is written to a C++, Python or Fortran file, language chosen by the user. 
For instance, the Fortran output would be: 

C This File is Automatically generated by ALOHA 
C The process calculated in this file is: 
C Gamma(3,2,l) 
C 

SUBROUTINE FFV1.0 (Fl , F2 , V3 , COUP , VERTEX) 

IMPLICIT NONE 

DOUBLE COMPLEX FlO) 

DOUBLE COMPLEX F2(*) 

DOUBLE COMPLEX V3(*) 

DOUBLE COMPLEX COUP 

DOUBLE COMPLEX VERTEX 

VERTEX = C0UP*( (F2(l)*( (Fl(3)*( (0, -l)*V3(l)+(0, 1)*V3(4))) 
$ +(F1(4)*( (0, l)*V3(2)+V3(3)))))+( (F2(2)*( (Fl(3)*( (0, 1) 
$ *V3(2)-V3(3)))+(F1(4)*( (0, -l)*V3(l)+(0, -1)*V3(4))))) 
$ +( (F2(3)*( (Fl(l)*( (0, -l)*V3(l)+(0, -1)*V3(4)))+(F1(2) 
$ *( (0, -1)*V3(2)-V3(3)))))+(F2(4)*( (Fl(l)*( (0, -1)*V3(2) 
$ +V3(3)))+(F1(2)*( (0, -l)*V3(l)+(0, 1)*V3(4)))))))) 

END 

Aloha's naming convention for generated routines is denoted by NAME_X, 
where NAME is the name provided by Ufo and X is either for an amplitude 
routine or the number corresponding to the off-shell particle for a wavefunction 



routine. The convention for the arguments is detailed in Appendix A 

To produce the necessary off-shell routines, Aloha loops over the list of 
particles and replaces the on-shell wavefunction with a propagator. For particle 
1, the incoming spinor is replaced by a propagator, computing^] 

Spinor(2) * Gamma(3,2, 1) * SpinorPropagator (1) * Vector(3) 

which is analogous to: 

FFV_1(^ 2 , e 3M , g, m 1} T 1 )=g* ^ ^ ■ % P + ™] e 3)1 . (2) 

p z — mf + %m\Y i 

FFV_1 is the wavefunction of the first spinor, and would be F0V in the origi- 
nal Helas language. The momentum for the propagator is computed by using 
conservation of momentum at the vertex. 



3 The propagator object is defined internally to Aloha. 
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Similarly, for the outgoing spinor, we have: 



FFV_2(V>i, e 3 ^ ma, T 2 ) = i — — 2 , • — F" * * 7 M ^2 e 3M . (3) 

p 2 — m| + ^m2^2 

The syntax used by Aloha reads: 

SpinorPropagator(2) * Gamma(3,2,l) * Spinor (1) * Vector (3) 

which can be related to the FIV routine in the original Helas language. 

Finally, Aloha also creates the routine associated to the off-shell vector. 
This routine returns a vector wavefunction from a pair of incoming and outgoing 
fermion wavefunctions, and it is written by ALOHA as: 

Spinor(2) * Gamma (3, 2, 1) Spinor (1) * VectorPropagator (3) 

analogous to the analytic expression FFV_3, which would be JIO in the original 
Helas language: 

vT - 

FFV_30/>i,7/> 2 ,#,m3,r 3 ) = fofypipi ■ -1-2 2 , - 3 r ■ (4) 

Following the original Helas convention, Aloha uses the unitary gauge for 
massive particles and the light cone gauge for massless particles. Notice that 
Aloha is not designed to create external particles routines, e.g., IXXXXX and 
OXXXXX for incoming and outgoing fermions, VXXXXX for vectors, SXXXXX for 
scalars and TXXXXX for spin-2 particles. There is a relatively small number of 
routines needed, and it is therefore reasonable to use a static librarjj^J 

In order to cover the full range of routines required for any type of interaction 
(and model), Aloha still needs two additional types of special routines: 

• Conjugate Routines: The presence of majorana fermions and of fermion 
flow violation, requires for some processes the presence of the conjugated 
version of the routines [21] not present in the Standard Model. In the case 
of the previous example, the Lorentz expression would have to be replaced 
by: 

C(-2,2) * Spinor(2) * Gamma(3,2,l) * 

Spinor(l) * Vector(3) * C(3,-3) 

In order to distinguish these conjugate routines from the standard ones, 
Aloha adds to their name the letter C followed by the number of the 



4 Thcse external particle routines are available in Python, C++ and Fortran languages in the 
Aloha package. 
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fermion pair that is conjugated. In this case the name of the conjugate 
routine would be FFVC1_X with X being either 0, 1, 2 or 3 depending on 
the number corresponding to the amplitude or off-shell particle. Note that 
_1 (_2) returns the wavefunctions associated to the incoming (outcoming) 
fermion. 

• Multiple Couplings: In the Ufo scheme each Lorentz structure is asso- 
ciated to a single coupling. In practice, this means that some interactions 
are linked to more than one Lorentz structure (for example, this is the case 
for the Z boson in the Standard Model). However, it is convenient for 
the readability and the compactness of the code to have one single routine 
per interaction. Therefore, Aloha defines a series of wrapper functions 
which call different routines associated by their couplings. The name of 
these wrapper routines are the concatenation of the name of the associated 
routines combined by an underscore. If all of the associated routines start 
with the same prefix (with letter S, F, V, T) then this prefix is only kept 
the first time. For example, the wrapper function of FFV1_X and FFV2_X 
will be called FFV1_2_X. 

In some cases, this method provides redundancy. Consider a four-Higgs- 
interaction. Due to the symmetry of the Lorentz structure, the four wavefunction 
routines will each have the same content. This type of symmetry is detected au- 
tomatically by Aloha and the computation is done only once. However, in order 
to simplify the work of the matrix element generator, all four routines are ex- 
plicitly defined and have their own names. The first one is computed and the 
remaining three are written as an alias of the first one. 

3.2. Resulting output 

The major advantage of Aloha comes from using it in combination with a 
Monte Carlo (MC) event generator. To this aim, Aloha offers two options. The 
first is to use Aloha as a standalone program. The user can call Aloha from 
the shell: 

./bin/aloha UFO.path [ -f Fortran I Python I CPP] [-o output.dir] 

where UFCLpath is the path to the UFO model directory, and the two possible 
arguments are 

• the language in which the routine should be written (default is Fortran). 

• the path for the output directory (which is by default UFCLpath/LANGUAGE, 
where LANGUAGE can be either Fortran, Python or C++). 

The second possibility is to internally link the Aloha code to a MC event 
generator. This option is only possible for a generator written in Python, and is 
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currently used by MadGraph version 5. For each process MadGraph requests 
the required set of helicity routines, which is produced by Aloha. This first 
implementation of Aloha inside a matrix element generator played an essential 
role in fully validating the package. 

4. Tests and validation 

The validity and precision of the original Helas routines has been known for 
many years. A natural choice for testing Aloha routines is therefore to compare 
them with the original Helas routines. Since the format of the Aloha is differ- 
ent from that of Helas, a line-by-line comparison would not work as a general 
validation test. Instead, validation is provided by the calculation of amplitudes of 
a particular process using the two sets of routines and comparing results. Below 
is a full description for the validation tests, as well as an example of the power 
and utility provided by such tests. 

4-1. The validation process 

MadGraph version 4 uses Helas routines to compute values for the square 
of the matrix element for a specific phase-space point of a known process. New 
routines can be validated by computing these values using both the original 
Helas and the new Aloha routines. 

The idea is to use two different versions of MadGraph to compute these val- 
ues. We employ MadGraph version 4.4.44, which uses the original Helas rou- 
tines, and compare the results with those from MadGraph version 5.1.1.0, which 
uses the new Aloha routines. 

The comparison is made by specifying a point of phase space (randomly gen- 
erated by RAMBO [25]) and a center-of-mass energy, after fixing properly the set of 
external parameters. The values for each matrix element squared are computed, 
and the difference is given by AD: 

a n - 2 v IM AL0HA )| 2 ~ l-M(HELAS)! 2 

| A4 (ALOHA) | 2 + |7W (HELAS) | 2 ' {) 

Due to effects of numerical precision, one can not expect exact agreement between 
the two codes. Instead, the test is considered to pass if values of the squared 
matrix element computed by Aloha and Helas agree to at least five orders of 
magnitude, i.e., AD < 1 x 1CT 5 . 

These tests are carried out for different processes based on the Standard Model 
(SM), the Minimal Supersymmetric Model (MSSM) and the Randall-Sundrum 
type I (RS-I) model, the later including spin-2 particles. Additional tests are 
made by checking gauge and Lorentz invariance of each process. We observe a 
great agreement for all processes, as illustrated in Tables [2] - [Tj and explained 
below. 
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4-2. Standard Model processes 

For the Standard Model, we test a very large sample of processes. Namely 
around 800 2 —7-2, 600 2 —7-3 and 50 2 — > 4 processes. These are chosen to 
ensure that all helicity routines have been covered. 

Some examples of these results can be found in Tables [2] and [3j In both tables, 
the first column shows the process chosen for the comparison. The second column 
shows the difference AD given by Eq. ^ for a specific point of phase space. 

4.3. MSSM processes 

One common approach to extend the Standard Model is the 4-dimensional 
Supersymmetric theory, based on the popular field description delineated by Wess 
and Zumino [26] . Within this category of Beyond Standard Model theories is the 
Minimal Supersymmetric Standard Model, which relates bosons to fermions by 
unifying symmetries in the most straightforward way [271 I2E1 [291 EO] ■ It addresses 
the Hierarchy Problem by adding new states whose contribution effectively re- 
move the quadratic divergences to the Higgs mass. It also provides a natural 
candidate for dark matter if the presence of an additional symmetry, R-parity, is 
assumed. On account of being a popular extension to the Standard Model, the 
MSSM has been extensively searched for at Hadron Colliders and several model 
implementations are available. We can therefore use these implementations to 
test Aloha within a theory which uses the full set of implemented subroutines. 
Testing Aloha for the MSSM is similar to testing it for the SM, since the Lorentz 
structures are for the most identical. However, the presence of majorana fermions 
and of fermion flow violations, require the presence conjugated versions of the rou- 
tines [23] not present in the Standard Model. The same validation tests used in 



Sec. 4.2 for the Standard Model are therefore performed for the MSSM. Approx- 
imately 1000 2 —7-2 processes and 7000 2 — > 3 processes have been tested. The 
results are conclusive and 100% agreement is obtained. A sample of processes 
checked for the MSSM is shown in Tables [4] and |U 

4-4- RS-I processes 

Several models with extra dimensions have recently been proposed. Examples 
of such theories include the Large Extra Dimensional model (or ADD) [3TI 132] 
and Randall-Sundrum (RS) models [331 133]- These models contain spin-2 particles 
allowed to propagate through higher dimensional space-time that couple to the 
SM particles in a very particular way. To complement the implementation of the 
ALOHA helicity routines, it is interesting to test and validate also spin-2 routines 

To this aim, a simplified version of the Randall-Sundrum type I model is used. 
Consider a warped 5-dimensional model in which the extra dimension is spatial 
and compactified in a line segment using the S1/Z2 symmetry. The propagation 
on the extra dimension is restricted to spin-2 particles, which will be perceived in 
4-dimensions as a tower of excited Kaluza-Klein modes of the same particle with 
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Table g SM 2 -> 2 results 


Table SM 2 -> 3 results 
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Table [I[ MSSM 2^2 results 
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increasing mass. We implement a reduced version of the RS-I model: it contains 
only the first mode of the spin-2 particle explicitly defined, which couples to SM 
particles. Validation tests, similar to the ones presented above for the SM and 
MSSM, are performed. Approximately 500 2 — > 2 and 2000 2 — > 3 processes are 
tested, and some examples can be found in Tables [6] and [7} Notice the first mode 
of the spin-2 particle is defined as y. The results are again in full agreement. 

4-5. Beyond the HELAS library 

An important achievement of Aloha is the ability to create an helicity routine 
for any type of interaction, even ones not found in the Standard Model. In this 
context, the Aloha package, via MadGraph 5, has already been used for some 
analyses. 

As an example we cite the Strongly Interacting Light Higgs model (SILH). 
Within this theory, the Higgs is the Goldstone Boson linked to a new strongly 
interacting sector. This model has been implemented into Ufo and presented 
in [11] . The authors provide a validation example where the partial decay width 
of the Higgs into two W bosons is computed. This decay is modified by a Non- 
standard Model interaction corresponding to the following Lorentz structure: 

VVS8 = Lorentz (name = 'VVS8' , 
spins = [3, 3, 1 ] , 

structure = 'P(1,3)*P(2, 1) + P(1,2)*P(2,3) 

- P(-l,l)*P(-l,3)*Metric(l,2) 

- P(-l,2)*P(-l,3)*Metric(l,2) ') . 

This structure corresponds to a dimension 6 operator (and should therefore be 
proportional to the vacuum expectation value). 

The value for the partial decay width is computed by MadGraph version 
5, which uses Aloha, and compared to results previously known. This example 
not only validates the Ufo model but also the routines automatically generated 
by Aloha. 

A second example can be seen in the recent study on effective operators that 
modify the top production and decay [35, 36J. These analyses are motivated by 
the precise result on the top sector provided by both the Tevatron and the LHC, 
and which allow one to constrain the higher dimensional operators. For instance, 
one of these operators is known as the chromo- magnetic operator [31)] : 

{HQ)a^THG% 
C = — 2 ^ + h.c, 

where A represents the cutoff of the effective theory. Adding such a term to 
the Lagrangian of the Standard Model leads to additional interactions shown by 
Fig. [2] These additional interactions are not covered by the traditional Helas li- 
brary. In collaboration with the authors of |35j, we check the Aloha results 
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diagram 1 NP=2, QCD=1 diagram 2 NP=2 




diagram 3 NP=2, 0ED=1 , 0CD=1 diagram 4 NP=2, QED=1 



Figure 2: Interactions induced by the chromomagnetic operator. 

against i) their analytical computation, and ii) their private implementation of 
the Helas routines inside the MadGraph version 4 framework. A very good 
agreement is again obtained. 

The implementation of Aloha inside MadGraph version 5, combined with 
the ability to automatically generate the Ufo model from a general Lagrangian 
via FeynRules, Aloha has been used and tested on a wide class of differ- 
ent models. The list of currently available models is presented on the Feyn- 
Rules web-page. These models include, for example, the Higss Effective Theory, 
the MUED model, the NMSSM and the RMSSM. 

4-6. Speed comparison 

In addition to the validity checks carried out to test the accuracy of the 
Aloha routines, one extra comparison is performed to check the speed of the 
new routines. Table [8] shows the time required to compute the squared matrix 
element using the original optimized Helas routines and these generated by 
Aloha. The third column shows the time ratio (Aloha/Helas) revealing that 
the automatically generated routines are typically within a factor of two in speed 
of the optimized routines. 

5. Summary and perspectives 

The automation of cross section calculations and event generation at leading 
order (LO), as well as at higher orders, from the Lagrangian level requires several 
non-trivial steps. In this paper, we have presented a code that automatically 
generates the helicity routines required by a MC event generator in several widely 
used formats (Fortran, C++ and Python). This is an important step in achieving 
a fully automatized chain from the Feynman rules, or the Lagrangian level, up 
to the event generation level. 

ALOHA represents a non-trivial extension of the current HELAS library which 
is limited to renormalizable interactions and fixed Lorentz structures. The code 
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Table ^[ Computational time (s) to produce 
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is extremely flexible and easily extendible. A first simple example is the imple- 
mentation of the complex-mass scheme [37] for the gauge-invariant computation 
of amplitudes with unstable intermediate particles (such as top, W and Z). As 
a further example, the current implementation of the Lorentz algebra assumes 4 
space-time dimensions, but it could be easily generalized to any other number of 
even dimensions, as the algebra is symbolic and its actual representation enters 
only at the final stage of the output writing. 

In addition, the implementation of the so-called R2 Feynman rules, which are 
needed for the automatic loop computation within the Cuttools method [38J, 
can be easily obtained. The necessary modifications for the routines to run 
on GPU's can be included, and the corresponding Heget [39, 40J library can 
be automatically generated for any BSM model. Finally, extensions to matrix 
element generator other than MadGraph, such as that used in Herwig++, is 
straightforward. Work in these directions is on-going. 
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argument 


symbol 


Python 


Fortran 


C++ 


mass 


M 


float 


real 


complex 


width 


W 


float 


real 


complex 


coupling 


Coup 


complex 


complex 


complex 


scalar wavefunctions 
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list (3) 


complex(3) 


complex [3] 


fermion wavefunctions 


F 


list (6) 


complex(6) 


complex [6] 


vector wavefunctions 


V 


list (6) 


complex(6) 


complex [6] 


tensor wavefunctions 


T 


list (18) 


complex(18) 


complex [18] 


Amplitude 


A 


complex 


complex 


complex 



All Fortran, C++ arguments are double precision. The arguments in 
Python are only informative, since Python is dynamically typed. How- 
ever, the list of wavefunctions should be initialized with the correct num- 
ber of elements. 



Poles Programme' Belgian Science Policy P6/11-P and by the IISN MadGraph 
convention 4.4511.10. 



Appendix A. Argument convention 

Since the Fortran, Python and C++ languages are distinct languages, there 
are small differences in the procedure used by each language to call the Aloha func- 
tions. Table A. 9 presents the different types of arguments, and Table A. 10 



presents various examples on how the routines should be called in each language 
and what is the equivalent analytical formula. 



Appendix B. Structure of the computation 

In order to perform the necessary analytical computations in Python, we 
design a dedicated analytical package. The fundamental object of this package is 
the Variable object, which contains three different pieces of information: 

prefactor The numerical value by which the variable is multiplied, [def ault=l] 

power The power of the variable, [def ault=l] 



name A string representing the variable. 
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Table A. 10 Syntax Examples 



VVSS.O 



analytical: 

Python 
Fortran/ C++ 



Coup * Vfr^VfSsS* 

A = WSS_0(V1, V2, S3, S4, Coup) 

WSS_0(V1, V2, S3, S4, Coup, A) 



FFV_0 



analytical: 

Python 
Fortran/ C++ 



Coup * F 2 T IM F 1 V» 

A = FFV_0(F1, F2, V3, Coup) 

FFV_0(F1, F2, V3, Coup, A) 



FFV_1 



analytical: 

Python 
Fortran/ C++ 



Coup * F 2 T^P(v£ 

Fl = FFV_1(F2, V3, Coup, M1,W1) 
FFV_1(F2, V3, Coup, Ml, Wl, Fl) 



FFVC1 



FFVCl.l 



analytical: 

Python 
Fortran/ C++ 



analytical: 

Python 
Fortran/ C++ 



Coup * FiCTl&FiVi 1 

A = FFVC1_0(F2, Fl, V3, Coup) 

FFVC1_0(F2, Fl, V3, Coup, A) 



Coup * wrjctp/Kf 

F2 = FFVC1_0( Fl, V3, Coup, M2, W2) 
FFVC1_0(F1, V3, Coup, M2, W2, F2) 



VVV1 2 



analytical: 

Python 
Fortran/ C++ 



V^Vi{Cou Pl Vl vp + Coup 2 * r% p ) 

A = VVV4_5_0( V1,V2, V3,Coupl, Coup2) 

VW4_5_0(V1,V2, V3, Coupl, Coup2, A) 



C stands for the conjugate matrix, T stands for the Lorentz structures 
associated to the vertex, Pf is the (fermion) propagator associated to the 
particle i. Other symbols are defined in Table A.9 Note in particular 
that the order of the arguments in the case of conjugate routines (FFVC1) 
follow the incoming/outgoing convention (and not the particle order). 
The last routine (VVV1_2_0) is an example of routines with multiple 
couplings, corresponding to two Lorentz structures. 
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For instance, 3x 2 is Variable (name =;, x' , prefactor=3, power=2). Variable 
objects with identical names are assumed to be referring to the same mathemat- 
ical variable. Therefore, addition and multiplication are performed as usual. For 
example: 

Variable (name= ' x ; , pref actor =3) + Variable (name= ' x ; ) 

= Variable (name= ' x ; , pref actor =4) 

In the case of a sum (product) between different variables, the sum (product) 
returns an object which is essentially the list of the summed (multiplied) variables. 

Variable (name= 'x ' ) + Variable (namely') 

= [Variable (name= ' x ; ) , Variable (namely')] 

With such a structure it is possible to represent any type of scalar expression. 

Non-scalar objects (like 7 M , ...) have basically the same structure, but 
slightly extended. There are two additional properties: lorentz_ind and spin_ind, 
which are the list of Lorentz and spinor indices linked to the object. Each object 
is related to a representation, i.e., each component should be linked to a scalar 
expression which can be a number, a variable or even an expression. 

Aloha inputs are based on Lorentz objects, and the calculations are per- 
formed in seven steps: 

1. Change of momentum convention in order to pass from the UFO into the 
HELAS convention. 

2. Multiplication of the Lorentz expression (from Ufo) to the wavefunction 
or propagator as required. 

3. Evaluation of the expression at the abstract level and simplification where 
appropriated (to adjust prefactors, etc. ). 

4. Use of the representation of each object in order to create a representation 
of the result. 

5. Simplification of each component. 

6. Factorization of each component to have an efficient way to evaluate the 
result. 

7. Writing the routine. 

Once the final expression is obtained, the result is passed to a specific routine 
which writes out the corresponding Helas routines. 

ALOHA primarily uses 4 files that correspond to the different steps of the 
computation: 

create_aloha.py: contains the driver of the code and the part needed to build 
the different routines associated with each lorentz structure. 

alohaJib.py: contains all of the code needed for the analytical computations 
and for the evaluation of the representation of the object. 
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aloha object. py: contains the definition of all basic objects. The current list 
of basic objects are displayed in Table [1} Adding a new object is straight- 
forward since it is enough to specify the index of the objects and the rep- 
resentation. 

aloha writers. py: contains the code that writes all Aloha functions in differ- 
ent languages (currently Python, Fortran and C++). 
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